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Abstract 

This is the fourth paper in our series of five in which we test the Master Constraint 
Programme for solving the Hamiltonian constraint in Loop Quantum Gravity. We now move 
on to free field theories with constraints, namely Maxwell theory and linearized gravity. Since 
the Master constraint involves squares of constraint operator valued distributions, one has to 
be very careful in doing that and we will see that the full flexibility of the Master Constraint 
Programme must be exploited in order to arrive at sensible results. 
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1 Introduction 

We continue our test of the Master Constraint Programme (Ij for Loop Quantum Gravity (LQG) 
ini 13 El which we started in the companion papers |21 El IH and will continue in [3] . The Master 
Constraint Programme is a new idea to improve on the current situation with the Hamiltonian 
constraint operator for LQG 0. In short, progress on the solution of the Hamiltonian constraint 
has been slow because of a technical reason: the Hamiltonian constraints themselves are not 
spatially diffeomorphism invariant. This means that one cannot first solve the spatial diffeo- 
morphism constraints and then the Hamiltonian constraints because the latter do not preserve 
the space of solutions to the spatial diffeomorphism constraint ^U]. On the other hand, the 
space of solutions to the spatial diffeomorphism constraint JO] is relatively easy to construct 
starting from the spatially diffeomorphism invariant representations on which LQG is based 
which are therefore very natural to use and, moreover, essentially unique. Therefore one would 
really like to keep these structures. The Master Constraint Programme removes that techni- 
cal obstacle by replacing the Hamiltonian constraints by a single Master Constraint which is a 
spatially diffeomorphism invariant integral of squares of the individual Hamiltonian constraints 
which encodes all the necessary information about the constraint surface and the associated 
invariants. See e.g. |2j for a full discussion of these issues. Notice that the idea of squaring 
constraints is not new, see e.g. |12|, however, our concrete implementation is new and also the 
Direct Integral Decomposition (DID) method for solving them, see |2j for all the details. 

The Master Constraint for four dimensional General Relativity will appear in |14j but before 
we test its semiclassical limit, e.g. using the methods of |15| I16j and try to solve it by DID 
methods we want to test the programme in the series of papers (21 El IH El ■ In the previous 
papers we focussed on finite dimensional systems of various degrees of complexity. This time 
we face free quantum fields squarely and will have to worry about the associated ultraviolet 
singularities. We study free Maxwell Theory and Linearized Gravity (formulated in connection 
variables) both of which are models with an infinite number of Abelean first class constraints 
linear in the momenta. 



2 Maxwell Theory 

The canonical formulation of Maxwell Theory on consists of an infinite - dimensional phase 
spaceA^ with canonically conjugate coordinates (Aa,E'^) and symplectic structure 

{Aa{x),A,iy)} = {E''ix),E\y)} = 0, {Eyx), Ahiy)} = eHt5{x, y) (2.1) 

where e is the electric charge and we are using units so that a = he^ is the Feinstrukturkonstante. 
In particular, the U{1) connection A has units of cm^^ while the electric field E has units of 
cm . 
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The clean mathematical description of Ai models on a Banach space £ j^. Here we 
will not need all the details of this framework and it suffices to specify the fall - off conditions 
of A, E at spatial infinity. Namely, in order that the canonical classical action principle be 
well-defined, that is, the action be functionally differentiable, the fields A, E respectively must 
fall off at spatial infinity at least as r~^,r~^ respectively. 

Consider now the infinite number of Maxwell - Gauss Constraints 



G(A) 



(fxK{daE'')=:<K,d-E>i 



(2.2) 



where A is a smooth test function of rapid decrease and f) = L2(Ill, d'^2;). Consider a positive 
definite operator C on independent of M. and define the associated Master Constraint by 



M := ^ <d-E,C-d-E>i^ 



(2.3) 



Obviously M = if and only \id-E = a.e., that is, if and only if G(A) = for all test functions 
A of rapid decrease. For the same reason for a twice differentiable function O on we have 
{{M, O}, 0}m=o = if and only if {d ■ E, 0}m=o = a-Sj hence if and only if {G(A), O}i\/i=o for 
all test functions of rapid decrease. 

Recall that the Maxwell Hamiltonian is given by (oc means equal to on the constraint 
surface) 

H = ^ j Sab {E^'E' + B'^B') (Xhj d^X 5ab ^ Pab (2-4) 
where the tranversal projector is given by 



{P± ■ f)a = fa- daA-'d'fb 



and 



-1/2 



E" 



(2.5) 



(2.6) 



Notice that we have used the Minkowski background metric in various places which distinguishes 
this model from the subsequent ones in [Jj. In the subsequent models we will treat background 
independent theories which excludes background dependent Fock space quantization methods 
which we will employ for this model. 

In order to solve the Master Constraint for Maxwell Theory on the usual Fock space it is 
mandatory to choose a nontrivial operator C in order that M becomes a densely defined operator 
defined via annihilation and creation operators. 

Let bj be any orthonormal basis of f) consisting of real valued, smooth functions of rapid 
decrease. The index set I in which the / takes values is countable and we could for instance 
choose the 67 to be Hermite functions. Next, consider the Hilbert space f}3 := with scalar 
product </,/'> [,3:= / d^x 5"-^ fafb- Let us define the functions 



b'l 


:= V- 






:= V- 


-A^'^dabl 


4? 


:= V- 


-A ^dabl 



(2.7) 



Notice the relations 



<bf\bf> 



< b 



f)3 = <b'i,b'j>i, 
>f)3 = <bi,bj >t,= 61J 



(2. 
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(3) 

Let us complete the bj to an orthonormal basis of f)3 by choosing some smooth, transversal 
and orthonormal system b^^\ b^p of rapid decrease, that is, d"'b^jj = 0, j = 1, 2. Now from (|2.8|) 

(3)/ 

it is clear that the longitudinal vectors bj do not form an orthonormal system, but they can 
be, as elements of f)3, expanded in terms of the orthonormal basis b^jb^ ,bP where only the 
expnasion coefficients for the b^ are non- vanishing. 

We thus find, using the completeness relation several times 

M = ^<d-E,C-d-E>tj 

= ^^<d-E,bi>i^ <bi,C -bjy^ <bj,d-E>^ 
I, J 

= I E < - 4'^' >f)3 <bi,C- bj >^ < bf, [z - z] >f,3 
i,J 

= f E E < b'S^bf >^3 <bi,C. bj >, < bf, bP >,,] X 

M,N I, J 

X <[z- z],b^p >f,3 < bP, [z - Z] >f,3 
= ^ E E < V^^^^bM, bi >(, <bi,C- bj >f, < bj, V^^^'^bN >f,] x 

M,N I, J 

x<[z- z],bP < bP, [z - z] >^3 
= <6J,^/^'/'c^/^'/V>f, <[^-^],6f >„3<&^^[^-^]>f,3 (2-9) 



If we would choose C = ^J—/S. then (|2.9() would simplify very much, however, this is impossible 
due to the boundary conditions on E which would imply that for this choice of C the integral 
(|2.l-i|) diverges logarithmically. However, notice that the operator Q on fj defined by 

Q = V^^^^CV^^^^ (2.10) 

is positive definite. 

We now come to the quantization of the system. We consider the kinematical operator 
algebra 21 generated from annihilation operators 

zP ■.= <bf,z>^, (2.11) 

and the corresponding creation operators given by their adjoints which are subject to the usual 
commutation relations 

i^'j, z'k] = (z'k)^] = 0, [z^j, (z'k)^] = aS^HjK (2.12) 

We represent this algebra on the usual kinematical Hilbert space TiKin of Maxwell theory given 
by the Fock space J- generated from the cyclic vacuum vector 17 defined by 

z^jn = (2.13) 

In terms of creation and annihilation operators, the Master Constraint Operator becomes 

M = ^ E Q^"" [^5 - i^m^ [4 - (4)^] (2.14) 

J,K 

Notice that 1)2. 13(1 is not normal ordered, hence it is a non-trivial question whether ()2.13|1 is even 
densely defined. 
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Theorem 2.1. 

The Master Constraint Operator 
operator. 



is densely defined if and only if Q is a trace class 



Proof of theorem 12.11 

Since the Fock space is the closure of the finite hnear span of finite excitations of the vacuum 
Q, M is densely defined if and only if M has finite norm. We compute 



MQI 



J,K,M,N 



a 



QjkQmn < a5jK - izj)\z^)^]n, a6MN - izlj)\z%)^]n > 

J,K,M,N 

jf QjKQMN[a^SjKSMN+ <^,z]^zj{zl^)\z%)^n> 



J,K,M,N 



"9-'^ E QjkQmnISjk^mn + ^jm^kn + Sjn^km] 



J,K,M,N 



a , 



[2Tr(Q2) + (Tr(Q))2 



(2.15) 



Both terms in the last line of ()2.15|) must be finite since Q is symmetric (even positive). The 
first term in the last line of 1)2. 15(1 is finite if Q is a Hilbert - Schmidt (or nuclear) operator, 
the second if Q is trace class. Since every trace class operator is nuclear, it is necessary and 
sufficient that Q he trace class. 

□ 



Now the algebra of trace class operators comprises an ideal within the compact operators. 
Compact operators are bounded, have pure point spectrum and every non-zero eigenvalue has 
finite multiplicity, hence zero is an accumulation point in their spectrum unless Q is a finite 
rank operator. A possible choice for C would therefore be for example 



C = e^^^¥^^-"' (2.16) 

where I is an arbitrary finite length scale. The exponential in 1)2. 16(1 is nothing else than minus a 
three dimensional harmonic oscillator operator, so its eigenvalues are := exp(— [| + ni +n2 + 
ns]), rij G No, hence its trace is explicitly e~'^/^(l — e~^)~^ and that of its square e~^(l — e~^)~^. 

Theorem 2.2. 

The choice \2. 1 6)) makes the integral \2.'J\) converge. 

Proof of theorem 12.21 
We expand 1)2. 3|1 in terms of coherent states ijjz for the Hamiltonian H defined by 

= g-|klP/2 e^"^"|0> (2.17) 

where 

Z"" = ^[X'' /I - ilpa], Za = ^[x^'/l + ld/dx''] (2.18) 

Using the overcompleteness relation for coherent states 

^ ^3 ^ ^z < ^z,- >i)= id(^ (2.19) 
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we find 



1 f cP'zd'^z f d?z'd?z' 



X 



^ 1 /■ d'^zd^z r d^z'd^z' 



O /ri ^^^^^ 

2 



X 

ll^lli /■ d^zd^z f d^z'd^z' 



2 Jc^ 7(^3 7r3 

-3/2 , ,1, r, , _H , ^ n M ^^-3/2 



x||^/^ 'd-V'.llfe [\<^z,e-''A'>i,\] t^-V'.'llfe 



II^IIL r d^zdH r d^z'd^z' 



2 J(e3 7([,3 7r3 

x^< V'z,\/^~V^ >f, [| < V'2,e"^V^' >f, |] -\/<Vvy^^~Vv^ (2.20) 

where in the third step we have used the Schwartz inequality. 

Since the classical electric field energy IL-^lia converges, it suffices to estimate the integrals 
in (|2.2U|) . Using the expansion of the coherent states into energy eigenfunctions it is easy to see 
that 

I < V^,e"^Vy >f) I = e~^^'^\exp{-^[zaZa + z'^z'^ - 2zaz'^/e])\ 
= e-^/'^exp{-^^\za- zX) exp{-^{l - ^)[zaZa + z'^zj) 

a 

< e-3/2exp(-l(l-i)[z,z, + zX]) (2.21) 
2 e 

Next, using the Fourier transform we get (up to a constant phase) 

^P^{k) = J d^xe-'^^'ip.ix) = (2/V^)3/2exp(-|H|2/V2 + ||fe||^P/2-i/ifc.z) (2.22) 
and have by using polar coordinates 



d^k ||fc||2;2_2;2fc.p 



||A:||(2vr)3 
Jo J -I 



/-3 1 





1 /"'IIpII dr _ . 

e 

TT 



< j-^ (2.23) 
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Thus we can finish our estimate by 



2 vr^ 



1-i [j d'xe-'^W-Wy [j 



v4M 

MiTilp —3/2 / 1 /'oo 



2 vr^ V e JO 

2 vr V e 

where the last integral resulted in a F— function. 

□ 

In what follows we will not need the specifics of Q, we choose any trace class operator such 
that (|2.9p converges. 

Having made sure that both the classcial and quantum Master Constraint are well-defined 
we can proceed to the solution of both the classical and the quantum problem. 

2.1 Physical Hilbert Space 

The full Fock space J- can be conveniently written in the form T = ® where = .F^'^^ 
contains the longitudinal excitations while = J^*-^^ ®T'^'^^ constains the transversal ones. The 
longitudinal Hilbert space in turn acquires the direct sum structure 



— ©r (3)^ ^, n\-\ (2.25) 

{"/ I/ex; E/"/ <oo {nK*>\ ^ ' 

where the overline denotes closure and the Hilbert spaces •^|^{3)| are one dimensional and can 
be identified with the multiples of the vector 

n ^ (2.26) 

For our purposes it is more convenient to write 1)2. 25() in a different way: Denote by S'(X) the set 
of all finite subsets J <ZT. For each J G ^(T) consider the subspace t\j of JFH consisting of the 
closure of the finite linear span of vectors such that nf^ > if / G ^7 and nf^ = otherwise. 
The Hilbert space J-^^ is the inductive limit of the family of Hilbert spaces J-j where directed 
set S{T) is equipped with the inclusion as partial order. 

Following the general programme of the direct integral decomposition (DID) we should now 
simply find a cyclic system 0„ of mutually orthogonal C°°— vectors for M so that the closures 
of the finite linear spans of vectors of the form M il„ are mutually orthogonal for different n. 
While we know that this can be done in principle, it turns out to be a hard problem and we were 
not able to find an explicit system O^. We leave it to the ambitious reader to come up with such 
an explicit solution for the present case. In what follows we therefore follow an indirect strategy 
which is also useful for other complicated field theoretic problems where it may be hard, if not 
impossible, to find an explicit solution {rin}- 

It will be convenient to choose the bj to be the eigenstates of the trace class operator Q with 
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positive eigenvalues Qj. For the example we discussed above the 6/ would be Hermite functions. 
Then for any J £ S{I) we define 

M^--(f)'EQH#^-(#^)^f (2.27) 

One can check that the family of operators H2.27() is not an inductive family with respect to 
the inductive limit structure introduced, basically because it does not annihilate the vacuum. 
However, the following holds. 

Theorem 2.3. 

The family of operators \2.27\ ) converges strongly to M on the dense set of Fock states. 
Proof of theorem 12. 3t 

To see this let i/j € J-^^ be a Fock state, that is, a state which is a linear combination of states of 
the form (|2.2(ij) . Then we find a minimal <ZT such that ip G T^j^ . Now consider any J with 
C J. Then ^ 

||(M-M^)^||^„ = <V^,(M-M^)2V;>^|| 
= < 17, (M - M >^ll 

= &{2 Y: + [ E Q^]'} (2-28) 



lei-J lei-J 



where in the second step we used that the operator M — Mj commutes with all the annihilation 
operators which create ^jJ from the vacuum 0, together with \ \ip\\ = 1. In the last line we used 
the previous computation (|2.15|) . 

Now since Q is trace class, given e > 0, tl^ €^ J^^^ we find Ji such that C Ji and such that 



(|2.28|) is smaller than e. Since the series in (|2.28|) is monotonously decreasing as ^7 ^ X, this 
holds also for all J\ d J which establishes the proof. We note that the label J\ depends on e 
but not on ^ itself but only on J^^ hence the limit is partly uniform. 

□ 



The idea to solve the constraint now rests on the following observation. 
Theorem 2.4. 

Let On be a sequence of self - adjoint operators with common dense domain D such that i) 
CLnip — > aV' foi" all tp £ D where a is another self ~ adjoint operator also defined densely on D. 
Suppose a) that x,y app{a). Then s — lim„_+oo En{{x,y)) = E{(x,y)). 

Proof of theorem 12.41 

Recall that a,i — > a in the strong resolvent sense provided that Rz{an) Fiz{a) strongly for any 
(and therefore all, see theorem VIII.19 of 17j z S C with ^(z) / 0. Here Rz{o) = (a — is 
the resolvent of a. By i) and ^7j> theorem viii.25a), o„ ^ a in the strong resolvent sense. By 
ii) and theorem viii.24b) the claim follows. 

□ 

The theorem applies in particular when 0^,0 have purely continuous spectrum in which case 
the convergence holds for all measurable sets. 

We can apply this theorem to our case because both assumptions i) and ii) are satisfied. 
Namely, by theorem 12.3)1 all the operators M j are densely defined on the dense set of finite 
linear combinations of Fock sates and converge there to M strongly. Furthermore, the operators 
Pj := i[z^^^ — [z^^^y] are mutually commuting, self - adjoint and have absolutely continuous 
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spectrum. Hence also the Mj-, M are mutually commuting and self adjoint (by the spectral 
theorem) and have absolutely continuous sectrum. It follows that s — limj^x Ej(B) = E{B) 
for any measurable set B where Ej, E denote the p.v.m. of M j- and M respectively. 

Let be any vector such that /Xf7„ is of maximal type (not to be confused with the Fock 
vacuum). For instance it may be the vector which actually defines the spectral measure ^ 
underlying a direct integral representation of subordinate to E. We then have 

— — ' = lim — ' 

x^o+ < nQ,E{{-oo,x])Qo > < wo(0),a;o(0) 

where the last equality holds /i— a.e. and ip,i{j',ujQ are the direct integral representations of 
^,^^',^0 respectively. 

Since X is countable it is in bijection with N and thus we take / = 1,2, .. for simplicity of 
notation. Consider any n, m and let ^'m^ ^Om be vectors which are finite linear combinations 
of Fock states with excitations at most up to label I = m and let En be the p.v.m. of M„ := 
(a/2)^E/=iP?- We calculate with /i[^„(5) :=< Oo,^„(S)Qo > etc. 

l^noix) n'^^ix) fif^^ix) f^Qo^ix) 

f^<f,^'ix)f^noix) - f^l,^'{x)t^no{x) ^ ^^l,^5"ix)^^noJx) - 

^^no{x)^i'^^{x) n^^{x)n'^^Jx) 
< ^, [E{x) - Enjx)]^' > finojx) + fi^,q,'{x) < n, [Enjx) - E{x)]no > 
Hn{x)[fino{x)+ < n, [En{x) - E{x)]no >] 

[m>ix) + f^'^,„^^no,noJ^) + /^no,no™-no(2^)+ < ^o, [En{x) - ^(x)]J^o >] 

^ (2.30) 



+ 



[nno{x)+ < no, [En{x) - E{x)]no >] 



Since the finite linear combinations of Fock states are dense, for given e > 0, we find 

some m and corresponding unit vectors ^m, ^m, f^Om such that ||^ — ^m\\, ||^' — ^mlli ll^o — 
^Omll < £• Since En{x), E{x) are projections with uniform operator norm at most unity we 
find, using the Schwarz inequality e.g. |/x^_^(^^iii'(x)| < ||^' — < e or \fJ'^_x^^^^ ^/{x)\ < 

\\^ - ^'mll < e. Next, for given x > 0, e > 0, ^,^',Qo there exists no = n{x,e,'^ ,Qo) such 
that \\[En{x) - E{x)]^\\, \\[En{x) - El.x)]'^']], \\[En{x) - E{x)]no\\ < e tor all n > uq. It follows 
e.g. I < ^, [En{x) — E{x)]^' > I < e. The absolute value of 1)2. 3Up can therefore be estimated by 

< ,{t'^oix) + \MMx)\ 



^lno{x)[^lQoix) - e] 
f^UoJx) + H^,^>Jx)\ 

il^Uoix) - e][/^no(a;) - 3e] 
< 2e 1 1 

/zn,)(x)-e fiQoix) /iQo(x)-3e 

^ T -^12 (2-31) 
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where we have assumed that, given x > 0, JIq we have 3e < Unoi^) ^-i^d used \fj.Qg{x)\, \fi\^,^^^l{x)\ < 
1. It is clear that given any 5 > we may choose e such that the last line of (|2.31l) can be made 
smaller than 6. 

We now wish to to calculate the approximation 

^^1^:^ (2.32) 

to (|2.29j) more explicitly. Let := 0^ (g) where {Q,'^)n is a cyclic system for M„ on the 
Hilbert space J-n which is the completion of the finite linear span of Fock states with excitations 
at most in the first n degrees of freedom and is a Fock basis of the orthogonal complement 
(^ri)"'"- We may identify 0^ with J7^q where cq = 1 is the vacuum of [J^n)-^. A given vector 
^ can then be written in the form ^ = ^ ^'^^(M„)r2^^ for certain measurable functions 
^ Na- 

Up to this point m, n are uncorrelated. Choose n > m. Then the measurable functions 
^mAfa' "^mNa'^OmNa Corresponding to ^'m, Qom respectively are in fact polynomials 

and moreover they vanish unless a = 0. Introducing a direct integral representation of J^' 
subordinate to M„ based on the formula (|2.32j) becomes 



lo d^i-{^)[ENPW\^OmNi^r] 

where p%{x) are the Radon - Nikodym derivatives of the spectral measures //^(i?) :=< f^^o' ^n{B)Q 
with respect to the total measure 

" = E„,„2-(««) = E„2-« 

where we exploited that E^i^B^ does not act on the 6q/. 

In order to compute ()2.33p in the limit x ^ we now will use a convenient choice of VL"^. 
These are simple modifications of the example of n mutually commuting constraints that we 
discussed in the companion paper [S] . The Hilbert space is unitarily equivalent to the Hilbert 
space ^2(5', d^c) where S is the set of real valued sequences {pi)YLi and /ic is a Gaussian measure 
with white noise covariance, that is /iG(exp(iFj)) = exp(— ^ Ylii //) where Ff = fiPi- This 
follows immediately from 

^icie'^f) :=< r?,e-S^^^[^i''-(^i'')^lj] > (2.35) 

and by using the Baker - Campbell - Hausdorff formula. Here, as before, denotes the Fock 
vacuum which in the p— representation is given by i^{p) = 1. We now choose the vectors 
according to the first example in 0. Introduce the variables yj := pi\/Qi, I = 1, ..,n so that 
~ i Yll=i a'^d let h'^i, I ^ Lg he the harmonic polynomials of degree s = 0, 1, .. in the 
variables yj with respect to the Laplacian in the yj, / = 1, .., n. Here Lg is a finite set of indices 
which depends on s. We define 

n^=(,,,)({w}) ■.= ni{{pi}) :=c:,/i,,({x/})exp(--5]y|)exp(-5]yfM) (2.36) 

1=1 1=1 

where is a normalization constant which is unimportant for what follows. As we showed in 
0, the M^Qf^ are dense in J^i Using that formally dficip) = H/li e"^"''^^"'' dpi/V^^, we 
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now compute the corresponding spectral measures 

J, df,G{p)o{x-M^)\n-i{p)\^ 



d^y e--2^"-^yi \h-,{y)\^ 



dt e{x - 1) t^+"/2-i 



T{s + n/2) Jj,^ 

= /^"o(^) (2.37) 

where in the second step we have integrated out all but the first n degrees of freedom, in the 
second we changed variables to yj thereby cancelling the Jacobean 13/=! ^/VOi ™ numerator 
and denominator, in the third we used that for given s there exists an SO{n) rotation which 
transforms hso into hgi and the rotation invariance of d'^y, M„ and after this the calculation is 
exactly the same as in 0- The total measure therefore becomes as in (3j 

, oo 

^"(^) = 2 ^ 2- V,o(x) = ^^^^^^ _ 0(x - t) e-*/2 (2.38) 

and the Radon - Nikodym derivatives were calculated to be 



^^'^ ^ (i/i"(x) r(n/2 + s) 7o 



(2.39) 



which close to x = behave as 

2r(n/2 - 1) 



r(n/2 + s) 



x'[{s + m/2-l) - x/2]Ki25sfi (2.40) 



Notice that Lq = {0} so there is only one harmonic polynomial of degree zero. Inserting this 
into (|2.33() . formula 1)2. 32(1 becomes 



which close to zero becomes (the measures ^([0, x]) cancel both in numerator and denominator) 



EsiPim'^nsii^) ^"W(0) ^ M/;^oo(0) M/"^mOO(0) 

E.//'^Ko)|f]^_,(o)|2 in^moom' ^ ■ ^ 

up to terms which vanish in the limit x ^ by the intermediate value theorem of Lebesgue inte- 
gral calculus (remember that the ^^gi etc. are polynomials and non - vanishing for finitely many 
s,l only, hence the integrand is actually a smooth function and dfj,^{x) = dxa^{x) for a smooth 
function (T„ is absolutely continuous with respect to Lebesgue measure as follows from (|2.39|) '). 
Thus, 1)2. 42() approximates 1)2. 29(1 as closely as we want, that is, for all x > 0,5 > 0,^',^'',ilo 
we find m > mo(x, (5, JIq) a-nd n > nQ{x,6,^,^' ,Qo) with n > m such that these two 
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quantities differ at most by 6. Since (|2.42|1 depends only on the coefficient ^'^oo('-*) 

m,n it follows that the physical Hilbert space in the approximation given by 6 in one to one 

correspondence with the one dimensional span of the vector fiJ^.^^Q ^^q, that is with C Since 

this holds for all 6, the physical Hilbert space selected by M from JFH coincides with C Since 
JT = ^11 it follows Hphys = as expected. 

Notice that this result is independent of the concrete choice of Qj, that is, independent of 
the choice of trace class operator. 



3 Linearized Gravity 

In this section we will consider the constraints of linearized gravity on Minkowski background. 
Because of the Minkowski background we can apply the same techniques as for the Maxwell- 
Gauss constraint. 

We will work with the (real) connection formulation of canonical gravity, i.e. the canonical 
pair of fields {Ai, -^Ej) on a three-dimensional manifold S, where Aa is an su(2)-connection 

and = y^det(g)e°' is a densitized triad for the spatial metric qab- We use a,b,c, . . . for spatial 
indices and k, . . . for su(2)-indizes. The latter ones are raised and lowered with 6^^ and 6ik, 
so we do not worry about the position of these indices, k = SirG Newton denotes the gravitational 
coupling constant and /3 the (real) Immirzi parameter. Canonical gravity in this formulation is 
a first class constraint system. The constraints of the full non-linear theory are the Gauss, the 
diffeomorphism and the scalar constraint: 

Gj = daEj + CjkiA^Ef 



C = r—r-A FL - (1 + f)ejn.nKTK^)e,uiEtE'i . (3.1) 

In the last line Ki is given by /? Ki = Aa — Ta, where Tq is the spin connection of the triad e". 
F = 2(dA + ^4 A ^) is the field strenght and the curvature of the connection A. 

We assume that the fields are asymptotically fiat and adopt the boundary conditions from 
dHl, that is Ai (and Ki) falls off as r'^ and E^ ~ ^^^E^ + 0(l/r) at infinity, where ^^^E'^ is a 
densitized triad for a flat metric. 

Linearized Gravity in the connection formulation was previously considered in jlOj , however 
there complex connection variables '^Ai = T^ + iKi were used. Since we have P Ki = Aa — Ta, 
we can express the complex connection in terms of the real variables as 

^A^^='-A^^ + {l-'-)Ti . (3.2) 
Moreover gave the ADM-energy functional in terms of the complex connection as 

H = - [ d^x 7dit(^(«IJ^Afe" - ^^Afe*) (3.3) 

where '^Aab = ^Aie^. This energy functional can be derived as the surface term one has to add to 
make the integrated Hamiltonian constraint C{N) = d^xNC (where N is the lapse function) 
differentiable. The differentiability is achieved with respect to the fields Ai,Ej fulfilling the 
asymptotic flat boundary conditions. 
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Rewriting the energy functional (|3.3|1 with help of into the real connection variables 
gives 

^ = ^ ^ d'xVd^[{l+(3')[T^T,^ - ry,'] + [A,'A,^ - a: A,'] - 2[T,'A,^ - T,^A,']] (3.4) 

where we introduced Aab = '^Aiel and Tab = ^ael- 

We come now to the linearization of the theory around Minkowski initial data Aa = and 
Ej = 6j. We will apply the following theorem, proved in a more general setting in j2Uj . which 
deals with the linearization of a finite dimensional first class constrained systems: 

Theorem 3.1. Let (M, fi) be a finite dimensional symplectic manifold with coisotropic con- 
straints Cct and Hamiltonian H. Suppose that niQ £ M is a point on the constraint surface at 
which the Hamiltonian vectorfield of H vanishes. Consider the tangent space V := Tmp(M) and 
introduce the linearized constraints, quadratic Hamiltonian, and linearized symplectic structure 
respectively by 



■■= {x",^'}Kn := {m^m^}(mo) = Q^^'imo) . (3.5) 

Here are understood as coordinates of a vector fields x in the linearized phase space V , given 
by the relation = x-m'^ where m"' are coordinates on the phase space M (seen as -functions 
on M). The statement is: 

The C^*" are coisotropic for the symplectic vector space (y, fi'*"), even Abelian, and 
is invariant under their Hamiltonian flow. Moreover, the linearization of the reduced theory 
coincides with the reduction of the linearized theory. 

We will use the prescription of the theorem for our infinte dimensional theory and check 
afterwards whether the statement holds also for this case. Thus we take mo = {Aa = 0, Ej = 6j) 
and consider the tangent space at this point, coordinatized by '™^a and , so that we can 
write Aa = ''"^a and Ej = + ^™-f/J. We insert the latter into (|3.1j) and keep only terms of 
linear order in {^™Ai,^™Hj). The result is 

= db''^Hba + eabc''''A,b (3.6) 
= dJ'''Abb-db''^Aab (3.7) 

C^™ = ^abcda^^^Afyc (3.8) 

where we introduced '^^a'' := ^^^AiSj and = and all indices are pulled with 

respect to the flat background metric 5ab- (Hence we will not worry about index positions.) The 
induced symplectic structure is {^^^Hab{x),^^'^A'^'^{y)}ii^ = [3K,5aSi\5{x,y). 

To compute the linearized Hamiltonian we notize that 1)3. 4() is at least quadratic in Tab and 
Aab- So we need Tab in terms of ('™Aa, '"-ffj) to linear order, which is 

linp^f, = ^e^'^^da^'-^H^d + dd{5a^^H,, - ''^Hac " ""^i/ca)) • (3.9) 

This can be computed by using the condition = daE^ + e^kiT^EI for the spin connection F^. 
We can then replace F, A, ■s/det{q) by ^™r, ^™j4, 1 respectively in order to arrive at the quadratic 
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Hamiltonian 



s 



Now it is merely a computational effort to check, that the linearized constraints ()3.6p are indeed 
Abelian and that the Hamiltonian (|3.1U|) is invariant (on the constraint hypersurface) with 
respect to the Hamiltonian flow of the linearized constraints. 

Summarizing, we use for the linearized gravitational field canonical variables {^^^Aab,^^^H°'^) 
(and drop in the following the superfix 'lin'), subject to the Abelian constraints ()3.6p . The 
boundary conditions for the linearized fields can be read off from the boundary conditions for 
the non-linear theory to be Aab ~ ''^"^ and H"''' ~ r~^. 

As in the Maxwell case we will use Fock space methods to quantize this system, so we have 
to introduce a pair of complex conjugated fields whose smeared forms will be finally promoted 
into annihilation and creation operators. Our choice is 



Zab 



V2f3L 



W-'/^ • {Aab + eacddcHdb) - im^l'' ■ H, 



ab 



and Zab , (3.11) 



where we introduced the abbreviation W = \/—A. The non-vanishing Poisson brackets are 

ih{Zab{x),Zcd{y)} = 6{x,y)dacSab- 

With this choice the Hamiltonian can be written as a manifestly positive function, as will 
be shown later. 
Now we plug 

Hab = ^W-'/^ ■ {Zab - Zab) =: IpW-'/^ ■ Pab 

Aab = ^W'/' . {Zab + Zab) " ^^-'^^ ' ^aMz^b " Zdb) 

=: f3lpW^'^ -Xab-lpW"^'^ -eacddcVdb (3.12) 
into the constraints 1)3. 6|) arriving at 



Ga := Ga-W~^ -daC -W-^ -dbeabcVc 



c 



pip W^/^ ■ {eacb - W-^ ■ edbcdadd + VF"' • ^adbdddc) ■ x^, + Ip W-^'^ ■ d^Sa 



zPbc 



Pbc 



Va = pip W'/^ ■ {daSbc - dc6ab) ' X^c + Ip W'^'^ ■ {eucdadd " euadcdd^ 

C = pij,W^/^-eabcdaXbc + lpW-^/^-{-W\,-dtdc)-Pbc (3.13) 

where we introduced a more convenient equivalent set of constraints. 

For the definition of the (classical) Master Constraint we use the Hilbert spaces t) = L2(R^, d^x) 
and [)^, with inner product on t)^ being given by < fa, fa >t)3-= f d^xd""^ fafb- Later on we will 
also need f)^ with < f'ah, fab >t)^-= f d?x6'^'^5^'^f'abfcd- (Here we deviate from the usual mathe- 
matical notation and write the indices in the inner product in order to keep track of the ordering 
of the indices. For f)'^ this is not necessary, therefore the indices are sometimes omitted.) Now 
we can define the Master Constraint as 

M = \{<G,Ki-G>^i + <V,K2-V >^^z + <C,K^-C>^) (3.14) 

where Ki,K2 are positive definite operators on f)^ and is a positive definite operator on f). 
These operators have to be chosen in such a way that the Master constraint is well defined (on 
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the classical phase space) and can be promoted into a densely defined operator on the quantum 
configuration space. In the following we will rewrite the Master constraint using the fields Xq}^ 
and Pah-, smeared with some basis of . To begin with the construction of this basis we introduce 
an orthonormal basis hi of \] consisting of real valued smooth functions of rapid decrease. From 
this we define a basis for the longitudinal modes in by b^f^ := ■ dabj and complete these 
to a full orthonormal basis of by choosing smooth, transversal functions b^^^ , bf'J of rapid 
decrease. Finally we equip f)^ with the following basis (the index i assumes values i = 1,2): 

41c •= (^acfc + W'"^ • eadbdddc) • 4a l^ft long. right transv. 

bfl := -W'^ ■ dbbf} = • dbdchi left and right long. 

4tc W^""*^ ■ ^c4b l^ft transv. right long. 

4?c := {dadbc - dc6ab)hf} = -2-1/2^-? {^6bc - dcdb)bi symm. transv., trace part 



,(7) ._ 2-1/2^-2 . eMcdaddbf^ = -2-'/^W-' ■ ebdcddbi antisymm. transv. 



(3.15) 



We complete this system to an orthonormal basis of f}^ by choosing a basis 6^^|,, bf^^ of the 
symmetric transverse traceless modes^. Furthermore we introduce another basis of the left 
longitudinal right transversal and left transversal right longitudinal modes: 

b^iic ■■= ■ db^acbfa = W-^ ■ e^daddb^iS. left long, right transv. 

4lc^ ■= • ^hdndcddbfa = -W-^ ■ ebdaddbf^l left transv. right long. . (3.16) 

Since the operator 5 : tc ^ ^cdadd^a squares to W"^ = —A on the subspace of the transversal 
modes in f)^, one can also write 

4;:^ = . e^daddb^:^ 

^ = ■ ^.^ddbfi^ . (3.17) 

If one chooses the basis of the transversal modes in if" such that 6^^^ = W~^ecdaddb^i} one obtains 

, (1'2) _ , (11) , (I'l) _ , (12) 

"ibc — ^Ibc "ibc — "ibc 

,(4'2) _ ,(41) _ J42) , . 

"ihc - -'^Ibc °Ibc - -"ibc ■ \'^-^^) 

Now we come to the calculation of the first term in the Master Constraint (|3.14() 

<G,K^.G>^S = Yl <Ga,4a>f,3<4a'^l-^Ja>D3<^Ja>G,>„3 .(3.19) 

I,J;i,j=l,2,3 



^The completeness of this basis is verified in appendix 1X1 
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Using the explicit form H3.13|l of the modified Gauss constriant we can write 

= Ip ( ^ < ^^fec, (eacfe + ■ eadbdddc) ■ h^jl >(,9 

where in the second line we used that i>f}^, b^^J are transversal and used an integration by parts. 
In the second last line we expanded the transversal vectors W^^'^ ■ z = 1,2 in terms of the 
basis b%j = 1,2. 
In a similar way 

<Ga,bfl>^, = lp<Pbc,W^^^-bfl>^9 

= lpY.<Pb^^^jL>^^<W'^'-bfa^bf^>^^ ■ (3.21) 
J 

Thus, utilizing 

I,J;i,j=l,2,3 

=<bf,lw'/'K,W'/'.b^j:i>^s (3.22) 

we get 

<Ga,Ki-Ga >^3 = ll J2 < ""be, >f,9 + < mc, b^n!l >f,0 < 4a ' W^'^'i^iH^'/' • 4'] >i,^ 

I,J;i,j=l,2 

(/? < 6|>Vi,X6, >^9 + < b^ji^lpbc >^9 ) 

+/2 ^ < >^s< bf},W'/'K,wy' . bf^ >,s< b%p,c >^9 (3.23) 

i,j 

where we assume, that Ki (and later also K2) commutes with the projector on the transversal 
modes. We will see that this is always possible to achieve. 

In an analogous way as for Ga we get for the Va contribution 

< K, K2 ■ Va >^3 = ll < 4S >f)« + < Pbc, 41? ) < bflw^/^K2W^/^ ■ b[jl >^3 

I,J;i,j=l,2 

(/? < b^jS,Xbc >f,9 + < b%{>,pbc >i,9 ) 

+21I (/? < >^9 + < ptc, b%l >^9 ) < bf^,W'/'K2W'/' . bf^ >^3 

I,J 

{P < b% Xbe >f,9 + < 4^i,P6e >f,9 ) (3.24) 
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To obtain the C contribution we need 

<CM>^ = lp<P W^/^ ■ eabcdaXbc + W-^/^ • ( - W^6bc - dbdc) ■ Pbc, hi > 

= -pip < Xbc, W^/^ • eabcdabi >f,9 +lp < Pbc, W-^/^ • i-W^dbc - dbdc) ■ hi >^9 

= -2^'%{I3 < Xbc, W^^Hfl >^9 + < Pbc, w'/Hfl >^9 ) 

= -2^2/^ ^ (/? < Xbc, >f)S + < Pbc, bfl >^9 ) < W^/\i, hi > (3.25) 
J 

arriving at 

<C,K3C>fj = ^<C,hi>^<hi,K2hj>^<hj,C 



I J 



2/2 ^ (/? < Xbc, 4?c >f)S + < Pbc, hfl >^9 ) < h, W^I^K^W''!'' ■ bj > 



i,J 



(/3 < b%l Xbc >(,9 + < hfl, Pbc >^9 ) (3.26) 

The expression of the constraints in terms of the fields Zab and Zab is quite comphcated, 
therefore we will perform a canonical transformation which will simplify the constraints. We 
will describe the canonical transformation in terms of the coordinates 

4"^ :=< 6^^^,X6c >f,9 and pf^ ■=<hf^^,pbc>[fi ■ (3.27) 

We define new coordinates x^f^ and by 



Xj 


= i=(/34-)+pf)) 


Pi = 


^/2V'^^ 


-xf^) 


^(3) . 

Xj 




~(3) 

Pi = 






Xj 


= -L(;34-)+pf)) 


Pi = 






Xj 


= ±iPxf^+pr) 


~(6) 

Pi = 


^/2V^' 




Xj 


=^^{px^p+pr) 


~(7) 

Pi = 






Am 




Pi = 


pf' . 




Xj 


= xf ) 





(3.28) 



Here, for instance x'f *^ stands for 



x'^p ■.=<b^^;p,xbc>^.= <b'jil^bc>^9<bf;:\b^jip^. , (3.29) 

and the sum over J,j reduces to just one term if one uses a basis with the property (|3.18p . In 
this case the canonical transformation restricts to finite dimensional subspaces, indexed by /. 
This gives new complex fields , defined by 

4") = i=(4-)-.4")) 4") = -L(4")+,pf)) . (3.30) 

The Master Constraint is then 



2IJ ^I Xj 



I,J;i,j=l,2 I,J 7,J;ij=l,2 

+2 Y QSf -f-f + 2 E ^3/. #4'0 (3-31) 
i,J i,J 
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where qJ*/], Q^^] and Q3/J are the matrix elements of Qi = W^I'^K^W'^I'^ , Q2 = W^I^K2W^/'^ 
and Qs = W^/'^K2W'^/'^ respectively. 

We come now back to our assertion that the ADM-energy functional can be written as a 
manifestly positive function in terms of the fields Zab and Zab- We will calculate the energy 
functional in a specific gauge and then generalize to the whole phase space. To this end we note 
that the constraint hypersurface is given by the vanishing of all the x'^f^ , I G X where a runs 
through all the values li, 3, 4i, 6, 7, that is all modes except the symmetric transverse traceless 
ones. An natural gauge condition is to require the vanishing of p'^\ I ^ I where the index a 
runs through the same values as above. Thus in this gauge the fields Aab and Hat are symmetric 
transversal traceless covariant tensors of second rank. (Hence this gauge is called STT gauge.) 
In the STT gauge formula (|3.9() simplifies to 

r'afe = -eacddcH'^b (3.32) 

where the prime indicates that the fields are in the STT gauge. Using that —eacddc squares to 
—A on transversal fields the energy functional (|3.1UI) can be written as 

= ^^^'^(/5'r;fcC + «6-C)«,-C)) 

= ^ ^ {(i^HU-m'ab + W + eacddA){A',b + eaefdeH'j,) 

= 2 j^d^x7abihV^-z'^h) . (3.33) 

Since the energy functional is a gauge-invariant function (on the constraint hypersurface) we 
can rewrite it in terms of the fields which are not necessary in the STT gauge by introducing 
the projector P^^-* onto the symmetric transversal traceless modes (see ()A.2|) ): 

H'- = 2 / J'.^.Un./^ . F<«) . + terms van.hmg „„ the constraint snrface. (3.34) 

The canonical transformation (|3.28j) leaves the STT modes unaffected, hence one can express 
the ADM-energy functional equally well in the Zab,Zab fields by simply replacing the old with 
the new fields. We conclude that the energy functional is a manifestly positive function. 

We will now quantize linearized gravity and examine the conditions on the operators Ki,i = 
1,2,3. As in the Maxwell case we quantize the theory by introducing the kinematical algebra 
generated from 

^t'^ = <b\tSbc>t,9 (3.35) 

and the corresponding adjoint operators (zj"^)^^. This algebra is represented on the Fock space 
^ generated from the cyclic vacuum vector Q. The commutation relations are the usual ones 

= [(II") )t, = (li'^))^] = 6u6-^ . (3.36) 

The Master Constraint Operator is given by 

^pfx^ n fe') _L 9 n (33) ~(3) ~(3) , ^ fe') ~(4i) 



— ^\ 2^I,J;i,j=l,2^1IJ -^I -^J "I" ^ 2^/, J Vi/J Xj Xj ^ 2^lj-ij=i^2^2IJ -^I 



J 



+'^Ei,jQSf ^f^xf + 2Zi,jQ3iji'^PiT] (3.37) 



where qJ'/], Q^^fj and Q31J are the matrix elements of Qi = W'^^'^KiW^/^ , Q2 = 1^3/2^2^^3/2 
and Qs = 11^3/2^^1^3/2 respectively. 
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A calculation completely analogous to (|2.15|) reveals that M is a densely defined operator on 
^ if and only if Qi, Q2 and Q3 are trace class operators on f)'^ respectively f). 
One possible choice for these operators is 

Qj = P^e-^P^ + pWe-^pW for j = l,2 

Q3 = (3.38) 

where h is the three dimensional harmonic oscillator operator and P-'-,pll are the projectors 
onto the transversal and longitudianl modes in iy' respectively.-fCi = p-^W~^/'^e~^W~^f'^P'^ + 
p\\W~'^/'^e~^W~^/'^P^^ and analogously K2 commute with the projector P-*-. 

The following calculation shows that the classical Master Constraint is well defined with this 
choice of the QiS. We consider only the part < Ga-,Ki ■ Ga > since the other terms can be 
treated in a similar way. As in the Maxwell case we will expand this term into coherent states 
ipz (|2.17|) for the three dimensional harmonic oscilator h. In the following calculation we will see 
the components (P-*" • G)a and (pH • G)a for a = 1, 2, 3 as scalar functions. To emphasise this, we 
will place a dot above the a: d. We will also replace the label z S (D^ in by Za S C^, d = 1, 2, 3. 
We can then write: 

< Ga, Ki-Ga>i,3 = E E < P'-"'^W-'^^Ga, P^-'^^W-^'^ G a >f, 

a=l,2,3 7=±,|| 

X < V.,,e-V.'. >f,< ^lJ,'.,P^^^W-^l^Ga >f, 

(fzad^Za f d^z'-d?z'- 



X < V..,e-'^^4 >f,< i>,'^,W-^d,dcP^"'^W-^/^Ga >(, 



a=l,2,3 7=±,|| 



X < V..,e-''V.' >f,< W-^d,i^,,,W-^d,P^'<^W-hGa >f, 



x||T^-3/29,V',J|^3||Ty-3/29,V'^,||^3| < ^.^e-'^V'^ >f, I (3.39) 
The first factor in the last line can be simplified to 

||^^-l/2g^p(7)vi/-l/2G',||j3 = \\{pi^)G)a\\l ^ WGaWl (3.40) 

since partial derivatives commute with the projectors P^'^\ The remaining integrals in ()3.39|) 
are the same as in the sixth line of 1)2. 2U() , so we can use the estimates ()2.2UMT^1)) , showing that 
the Master constraint is a well defined phase space function. 

If we assume,as is the case for the Qj's in (|3.38() . that it is possible to choose 6/ and hf^, b^^^ 
such that 67, b^j]^, bf'J and b^^^ = ■ dabj are eigenstates of Qi, Q2 and Q3, we can write the 
Master Constraint Operator more compactly as 

M = |E'^i"^(^f^)' ' (3-41) 
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where the index a assumes values a = li, 3, 4i, 6, 7 with i = 1,2 and the qf are positive eigen- 
values (multiplied by two for a = 3, 6, 7) of the Q^-'s. Thus this Master Constraint Operator has 
the same structure as the Master Constraint Operator for the Maxwell theory H2.14() . Going 
through the same procedure as in this case, see section 12.11 we see that the physical Hilbert 
space is unitarily equivalent to the Fock space generated from the vacuum Cl by applying just 
the creation operators generating symmetric-transverse-traceless modes. 

4 Conclusions 

The basic lesson that we have learnt from the examples studied in the present paper is that, 
surprisingly, the Master Constraint Programme can cope with the UV singularities that one 
expects from squaring operator valued distributions. Take the Maxwell case for simplicity. The 
infinite number of classical Gauss constraints G{x) = daE'^[x) can classically be encoded in 
the single Master contraint M' = ^ J d^xG{x)'^. This is classically well defined because G(x)^ 
decays at infinity as and, moreover, it has well defined distributional second derivatives with 
respect to the phase space differentiable structure. However, quantum mechanically, M = oo is 
hopelessly divergent. This cannot be even repaired by subtracting an infinite constant because 
M , in contrast to the electromagnetic field energy H = ^ J d^x5ab5ab{E°'{x)E^{x)+B"-{x)B^{x)) 
has a different structure in terms of creation and annihilation operators. This is best seen by 
using the Fourier transform of the annihilation operators Za{k) = J S'xl\f2/n e~^^'^ Za{x) where 
Za{x) = ix/^Aaix) - z( v^)"^^"(x))/ ^2^. Then we get 

H = nj d'k\\k\\{6'''-^)[za{k)z,{k)^ + za{k)^z,{k)] (4.1) 
m' = I y d^k \\k\\ k'^k' [za{k)zbik)^ + Uk)^h{k) - Za{k)Zb{-k) - Za{k)Hb{-k)^] 

Thus, the normal ordered expression : H : is densely defined and positive while : M : is neither 
densely defined nor (formally) positive. 

We used the flexibility in the definition of the Master Constraint in order to circumvent 
this problem. What we did is to introduce a positive integral kernel K(x, x') (that is, a 
positive operator on the one particle Hilbert space \) = L2(K^,d^x)) and to define M = 
J d^x j d^x' K{x, x')G{x)G{x') which still classically encodes all the constraints. The idea is 
that the kernel serves to smoothen the operator valued distributions into well behaved opera- 
tors. We found that the necessary and sufficient condition for M to be densely defined is that 
K be trace class. For any such choice of K the physical Hilbert space selected by the Master 
Constraint from the full Fock space of all modes was the correct one, the Fock space for the 
transversal modes. Thus, the physical Hilbert space constructed is independent of the concrete 
choice of K which matches nicely with the fact that the classical constraint surface defined by 
M = is independent of the choice of K. 
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A Completeness of the basis 

Here we will verify that (|3.15|) is indeed a basis of the Hilbert space of square integrable covariant 
tensors of second rank t)^ = {L2(R'^))^. We will begin by establishing, that the modes listed 
there are complete. To this end we introduce the projector p onto the transversal modes in 

(p ■ v)a := pi ■ Vb := {6i + W-^ ■ dadh) ■ Vh (A.l) 
and define the following projectors on [)® 



(p{l) 


■TU 


= - fa) 


■pi- 


Ted 




2 left long, right transv. modes 


(p{3) 


■T)ab 


= {^l - Pi) 


■{St 


-pi) 


■Ted 


1 left and right long, mode 


(p{4) 


■T)ab 


=pl-{st- 


Pt)^ 


Ted 




2 left transv. right long, modes 


(p{6) 


■T)ab 


= \pab-p"' 


■ Ted 






1 symm. transv. trace part mode 


(p{7) 


■TU 


= Wa-pi 


-Pi 


pi)^ 


Ted 


1 antisymm. transv. mode 


(p{8) 


■T)ab 


= Wa-pi- 


^pI 


pi- 


PabP"") ■ Ted 


2 symm. transv. tracefree modes 



Using the projector property p ■ p = p, it is easy to see that the projectors are orthogonal 
to each other and satisfy • p(^) = (5"/^p(") . For instance 

(p(6).p(8).5^)^^ = hafP'^^UPl-p'f+P}-pi-PefP"')-Tcd 

= iPa6(/"+/"-/^<Je//")=0 (A.3) 

where we used p^^ 5ef = 2. Furthermore the sum of all the projectors in HA.2|) is the identity on 

Next we have to show that the h^f^ in (|3.15() are a complete basis of the subspace of a-modes 
in . We begin with the left longitudinal and right transversal modes, rewriting them as follows 

41? := {^acb + W-^ ■tadbddde)^hf^ = -W-^ -dbeedaddbf, . (A.4) 

This equation can be proved by using that every transversal vector tia can be written as the 

curl operator applied to another transversal vector tia = ^abcdbt2c- Additionally one utilizes 

that eabcdb ■ (-cdedd = —A ^ae on the subspace of transversal vectors. Now we have for any left 
longitudinal tensor Tab 

Tab = {SI - Pi) ■ Teb = W-'da ■ (-H^"^ • deTeb) =: W'^ ■ Oafb (A.5) 

SO it can be written as a gradient of a vector VF"^ • /{,. If Tab is additionally right transversal 
this vector has to be transversal. Therefore {bjed = • datdj}je:3 is a(n) (orthonormal) basis 
for the left longitudinal right transversal modes if {tdjjjex is a(n) (orthonormal) basis for the 
transversal modes in f)^. Now if the latter is the case then {t'^j := ■ eabddbtdj} j<^z ^^^o 
an orthonormal basis for the transversal vector modes, since the operator ta ^ ■ eabddbta is 
unitary on this subspace (on the transversal modes it squares to the identity and on the whole 
space to the projector p). So we have proved, that 6^^*^ and b^"^'^^ is a basis for the left 
longitudinal right transversal modes. In a similar way one can deal with a = 3,4i,4'i. 
For the completeness proof of 

= ■ {da6te - deSab)bf^ = 2-1/2(5,, + ■ dedb)hf (A.6) 

b\^l=W-''-aabi 

we write for a tensor Tab G P^^\^^): 

Tab = {P^^^ ■ T)ab = kPabP^'^Ted = \ {Kb + W^"' " Qadb) ■ Tee (A.7) 
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where in the second equation we used that T^t is right and left transversah T^c is a square 
integrable function, so it can be expanded into the basis of I}. We conclude, that 

i^7bc}-fe3 is a complete basis of the symmetric and transversal modes with tracepart. 

To verify the completeness of the antisymmetric and transversal basis b^j^^ = —2~^/'^W~^ ■ 
^bdcddi>i, we use that every antisymmetric tensor Tab can be expressed as 

Tab = f^abctc where tc = \ecdeTde ■ (A.8) 

The transversality condition daTab = translates into the vanishing of the curl eabcdatc of to so 
that tc has to be longitudinal, i.e. can be written as the gradient of a function ■ f. If Tab is 
square integrable, so is tc and hence /, so / can be expanded into a basis 6/ of f). We conclude, 

(7) 

that {bjf^c}i<^3 is ^ complete basis of the antisymmetric and transversal modes. 
This finishes our verification, that the basis 1)3. 15(1 is complete. 
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